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INTRODUCTION 


The  diffraction  model  allows  for  a  consistent  treatment  of  diverse  scat¬ 
tering  phenomena.  For  example,  when  electron-phonon  interaction  is  the  scat¬ 
tering  mechanism,  the  diffraction  model  provides  a  unified  treatment  of 
electrical  and  thermal  transport;  phonon  contributions  to  the  electron  effective 
mass;  and  acoustic  attenuation  in  perfect  and  point-  or  line-defected 
crystalline  alloys,  and  also  in  liquid,  amorphous,  and  disordered  alloys. 

The  essence  of  diffraction  model  methods  is  the  factorization  of  scattering 
rate  matrix  elements  into  terms  such  that  the  lattice  vibrational  properties  are 
entirely  contained  in  dynamical  structure  factors,  which  can  be  related  to 
transforms  of  generalized  pair  distribution  functions. 

The  most  celebrated  dynamical  structure  factor  is  tnat  defined  by  Van 
Hove  (ref  1) ,  viz. , 

S(K,fl)  =  /  dt  exo(iflt)  <  r(K,t)r+(K,0)  >-j  (i) 

-00 

where  <0>T  denotes  thermal  average  in  the  equilibrium  ensemble  (i.e.,  <0>j  = 

I  <i|  0  |  i>),  and  r(K,t)  is  the  spatial  Fourier  transform  of  the  ion  density 
operator.  The  Van  Hove  dynamical  structure  factor  S(K,ft)  describes  very  general 
scattering  phenomena  in  which  the  lattice  absorbs  pseudo-momentum  K  and  energy 
ft.  It  is  directly  measurable  via  neutron  scattering  experiments.  It  encom¬ 
passes  Debye-Waller  factors,  umklapp  processes,  couplings  to  transverse  phonons, 
multiphonon  effects,  etc.  Thus,  S ( K , S2 )  has  been  the  object  of  intense  theoreti¬ 
cal  interest. 

Baym  (ref  2)  elaborated  a  diffraction  model  approach  to  the  treatment  of 
electron-phonon  interaction  controlled  phenomena  based  upon  the  Van  Hove  dynami¬ 
cal  structure  factor.  Baym's  approach  has  proven  to  be  extremely  productive  and 
has  become  standard.  However,  Baym's  form  of  the  diffraction  model,  based  on 
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the  Van  Hove  dynamical  structure  factor,  is  strictly  correct  only  for  the  scat¬ 
tering  of  plane  (or  pseudoplane)  waves  (PWs).  To  apply  diffraction  model  tech¬ 
niques  to  electron-phonon  interaction  controlled  processes  in  general,  it  is 
necessary  to  derive  diffraction  model  expressions  appropriate  for  treating  the 
scattering  of  Bloch  waves  (BWs). 

Therefore,  the  principal  objective  of  this  report  is  to  present  a 
coherent  exposition  of  a  generalized  diffraction  model  which  applies  to  PWs  or 
BWs  in  perfect,  disordered,  or  defected  crystalline  alloys  and  to  PWs  in  liquid 
or  amorphous  alloys.  This  aim  is  realized  through  the  introduction  of  a  gener¬ 
alized  diffraction  model  which  applies  to  BW  electron-phonon  interaction 
controlled  processes  and  reduces  to  the  standard  diffraction  model  ’n  pw  scat¬ 
tering  cases. 

An  important  secondary  objective  is  to  reveal  the  diffraction  model  origins 
of  the  approximate  expressions  appearing  in  the  standard  literature.  In  par¬ 
ticular,  the  standard  equations  used  to  describe  electron-phonon  interaction 
controlled  processes  are  approximate  diffraction  model  expressions. 

THE  GENERALIZED  DIFFRACTION  MODEL 

The  development  of  diffraction  model  factorizations  for  scattering  of  BWs 
is  demonstrated  in  this  section.  The  essential  diffraction  model  ideas  are  con¬ 
tained  in  Reference  1.  The  theoretical  development  follows  that  presented  in 
Reference  3.  The  methods  employed  here  were  also  used  by  Van  Hove  (ref  1).  The 
crucial  step  in  either  reference  is  the  application  of  Bloch's  Theorem  (ref  4). 

In  order  to  streamline  the  notation:  (1)  Let  q  stand  for  (q,j)  where  q  is  a 
phonon  wave  vector  and  j  is  a  phonon  branch  and  polarization  index.  (2;  Let  k 
stand  for  (k,s)  where  k  is  the  scattered  particle  wave  vector  and  s  is  an 
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appropriate  index.  For  example,  for  electron  states,  s  might  represent  the 
spin.  (3)  Consider  pure  materials.  (The  generalization  required  to  treat 
alloys,  via  the  introduction  o *  partial  structure  factors,  is  straightforward 
and  well  known.)  (4)  The  following  shorthand  notations  are  adopted: 

K  =  k-  k',A  =  KJ-G-H,  and  A’  =  K  +  G’  -  H’ 
where  G,  H,  G',  and  H1  are  reciprocal  lattice  vectors  (RLVs). 

Scattering  from  an  initial  state  containing  a  BW  of  wave  number  index  k, 

|  k,i>  =  ^  a(G)  exp[i (k+G) .x]  i  i>  (2a) 

G 

into  a  final  state  containing  a  8W  of  wave  number  index  k’, 

j  k',f>  =  Y  £(G)  expf i (k'+G) .x]  |  f>  (2b) 

G 

is  treated.  The  sums  are  on  RLV  G. 

The  perturbation  H'  producing  the  scattering  is  assumed  to  be  of  the  form 

H'  =  Y  [V(x-xt)  -  U(x-i)]  (3) 

1 

where  X]  is  the  position  of  the  ion,  which,  in  the  absence  of  the  lattice  vibra¬ 
tion,  would  be  at  position  1,  V(x)  is  the  ionic  potential,  and  U(x)  is  the 
unperturbed  potential.  We  think  of  these  potentials  as  fully  screened,  etc., 
and  use  the  same  letters  to  denote  the  Fourier  transforms, 

{ U ( A ) ,V(A) }  =  /  dx  exp[iA.x] {U(x) ,V(x) }  (4a) 

The  origin  is  chosen  in  such  a  way  that 

A.  1  =  K.l  +  (integer)27r 

We  also  denote  the  difference  in  the  transformed  potentials  (i.e.,  the  transform 
of  the  undisplaced  potentials)  by 
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W(A)  =  V(A)  -  U(A) 


(4b) 


Thus,  using  Eqs,  (2)  to  (4), 

<k',f|  H'l  k,i>  =  ^  <3(G)*a(H){V(A)<f|  exp[iA.x-|]|  i>  -U(A)exp[iK. l]<fj  i>)  (5) 

1,G,H 


and 

<k'  ,fj  H'j  k,  i>exp[i  (e-j-ef )]  =  /  j3(G)*ct(H)  {V(A)<fj  exp[iA»x-|  (t)]j  i> 

1,G,H 


-U(A)exp[iK. l]<f|  i>{  (6a) 

where  x-|(t)  is  in  Heisenberg  representation  and  we  have  used 

Hol  ^  -  eii  (6b) 

and 

Hol  f>  =  £f|  f>  (6c) 

It  then  follows  that  the  scattering  rate,  neglecting  factors  accounting  for 
occupation  of  the  states  k,  k1,  etc.,  is  given  by 


X(k,k',Q)  =  ^  Pi  |  <k’,f|  H’|  k ,  i >|  ^  6(6i-ef-Q) 
i  ,f 


=  I  I  il  Pi  I  <K' -fj  H'J  k ,  i >j  2  expti  (e-j-ef-Q) t] 
i ,  f 

=  ]>  /J(G)*a(H)[/S(G,)*a(H')]* 

G,H,G' , H 1 

{V(A)V(A')*S,(A,A,.fl)  +  W(A)W(A' )*a(K)6(ft)}  (7) 

where  the  geometrical  structure  factor  a(K)  is  defined  in  the  usual  way,  i.e., 

a(K)  =  ^  exp[iK. (1-1 ' )]  =  N  ^  exp[iK.d]  (8) 

1,T  d 
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The  dynamical  structure  factor  S’ (A, A', ft)  is  defined  as 

S’ (A,  A' ,ft)  =  S(A,A*,fi)  -  a(K )S(fl) 


(9a) 


where  the  dynamical  structure  factor  S(A,A',ft)  is  defined  as 

S(A,A',Q)  =  ^  expfiK. (1-1 ’ )] 

1,1! 


j  “  exp[-iflt]  <exp[-iA’ .U]  ’  (0)]exp[iA.u-|  (t)]>j  (9b) 

and  u -j  ( t )  is  the  displacement  of  the  ion  at  1  due  to  the  lattice  vibrations, 
i.e.,  x-j(t)  =  1  +  u  i  ( t ) .  Note  that  Eq.  (7)  allows  for  arbitrary  combinations  of 
disorder  and  defect  scattering;  for  example,  Eq.  (7)  allows  for  straightforward 
application  to  scattering  in  impure  and/or  point  defect  containing  crystalline 
solids  at  finite  temperature. 

S(A,A',S2)  is  a  generalization  of  the  Van  Hove  dynamical  structure  factor 
(ref  1),  i.e., 

S(K,K,G)  =  S(K,ft)  (9c) 

ard  S' (A, A' ,fl)  is  a  generalization  of  the  modified  Van  Hove  dynamical  structure 
factor,  defined  by  Baym  (ref  2)^  i.e., 

S' (K,K,8)  =  S' (K,n)  (9d) 

Generalizing  the  development  in  Reference  3,  one  obtains 

<exp[-iA' .U] ' (0)]exp[iA.U](t)]>y  =  exp[A' .2(1-1 ’ ,t) .A]  (10a) 

where 

Z(d,t)  =  -  ^  Wq  (<2nq+l>-<nq+l>exp[i9cjf  t]-<nq>exp[-i0cjr|:] }  (10b) 

q 


where 

wq  ~  eqeq/2NMQq  (10c) 
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and  d  =  1-1 '  is  an  ion  spacing  in  the  absence  of  thermal  vibrations;  N  is  the 
number  of  ions;  M  the  ion  mass;  Slq,  nq,  and  eq  are  the  phonon  frequency,  number 
density,  and  polarization  vector  at  q,  respectively,  and 

°d,t  =  q*d  +  ^qt  (10d) 

Note  that  Z(d,t)  is  a  symmetric  second  rank  tensor  which  is  independent  of  K,  A, 
and  A’ . 

Thus,  we  can  express  the  dynamical  structure  factor  S(A,A',fl)  in  the  form 

S(A, A* ,«)  =  N  ^  exp[iK.d]  J  ~  exp[-iQt]exp[A. Z(d, t ) . A ' ]  (11a) 

d 

=  ^  [A..  .A]..  . In(K,ft) . .  .[A'..  . A ’ ] /n !  (lib) 

n=0 

where  [A..  .A]  is  an  n-fold  open  product  of  the  vector  A.  Equation  (lib) 
defines  the  structure  tensor  In(K,fi)  of  rank  2n  and  is  obtained  by  expanding  the 
exponential  in  Eq.  (11a)  and  performing  the  sum  over  d  and  the  t  integration. 
(N.B.,  the  tensor  In(K,Q)  is  not  uniquely  defined  by  Eq.  (lib);  however,  the 
most  symmetric  forms  with  respect  to  the  left  and  right  inner  products  are  the 
obvious  choice  and  are  intended.) 

Combining  Eqs.  (7)  and  (11), 

X(k,k',fi)  =  ^  |3(G)*a(H)  [/3 ( G 1  )*a(H'  )  ]* 

G,H,G' ,H' 

(V(A)V(A' )*S' (A, A’ ,fl)  +  W(A)W(A')*a(k)6(Q)} 

=  2  vn(k,k')*..  .In(K,fl)..  ,Vn(k,k')/n!  +|w0(k,k')|*  I0(K,Q)  (12) 

n=l 
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where  we  have  defined  the  scalar 


W0 ( k , k ' )  =  <k '}  W ( x ) |  k>  (13a) 

and  for  n=l  to  »,  the  rank  n  tensor 

Vn(k,k' )  =  <k'j  ( [  ( iV) . .  .  ( iV)  ]V(x)  )|  k>  f  13b) 

and  we  have  used  the  relation 

^  j3(G)*a(H)  [A...AJV(A)  =  f  dx{^  j5(H)exp[i  (k'+H)  .x]  }* 

G,  H  H 

r(-V)..  .(iV)]V(x)>{^  £(G)exp[i(k'+G).x]} 

‘J 

=  <k'!  ( [  ( i  V) . .  .(-5V)]V(x)>  ]  k>  =  Vn(:<)  (13c) 

All  the  dynamical  properties  of  the  lattice  are  contained  in  the  structure  ten¬ 
sors  (In(K,fl)  !  n=0,  »} . 

Equation  (12),  with  (In(K,o)  j  n=0,  »)  defined  in  Eas.  (10)  and  (11), 
W0(k,k')  defined  in  Eq.  (13a),  and  the  { Vn ( k , k ' )  |  n=l,co)  defined  in  Eq.  (13b), 
constitute  the  generalized  diffraction  model  expression  and  are  the  principal 
results  of  this  report.  N.B.,  to  obtain  actual  scattering  rates,  one  must 
include  appropriate  factors  (e.g.,  (2rr/h)f (k) [1-f (k' ) ] )  in  Eq.  (12). 


PLANE  WAVE  SCATTERING,  n-PHONON  STRUCTURE  TENSOR  $n(K,fl),  DEBYE-WALLER 
TENSOR  W,  AND  THE  VAN  HOVE  DYNAMICAL  STRUCTURE  FACTOR  S(K,fl) 

Since  the  first  term  in  Ea.  (10b)  is  independent  of  d  and  t,  it  is  useful 

to  define  the  tensor 


Z'(d,t) 


^  wq  [<nq+l>exp(i9d/t)+<na>exp(-i0d(t)] 

q 


(14a) 


and 


Z 


Wr 


<2nq+l> 


q 


(14b) 
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Thus,  Eq.  (10b)  can  be  rewritten  in  the  form 

Z(d,t)  =  -  W  +  V  (d,t)  (14c) 

and  since  the  Debye-Waller  exponent  2W(K)  is  simply  related  to  W,  viz., 

2W(K)  =  K.W.K  (14d) 


we  designate  W  as  the  'second  rank)  Debye-Waller  tensor. 

Equations  (14a)  through  (14d)  suggest  that  a  different  factorization  of  the 
scattering  rate  will  be  particularly  appropriate  in  the  PW  case.  The  treatment 
of  PW  scattering  leads  in  a  natural  way  to  the  introduction  of  n-phonon  struc¬ 
ture  tensors  (4>n(K,ft) > ,  which  are  also  useful  in  expressing  limiting  BW 
expressions.  Thus,  Eq.  (2a)  becomes 

|  k,i>  =  exp[ik.x]  |  i>  (15a) 

etc . ,  and 


A,A'  — >  K 

Thus,  recalling  Eq.  (9),  Eq.  (11)  yields 


(15b) 


S(K,fi)  =  N  ^  exp[iK.d]  j  ~  exp[-iJ2t]exp[K.Z(d,t)  .K]  (16a) 

d 

=  ^  (K...KJ...  . In (K,fi) . . .  . [K. . .K]/n!  (16b) 

n=0 

Employing  Eqs.  (14c)  and  (14d),  Eq.  (16a)  yields 

S ( K , fi ) =exp [ -2W ( K ) ]  N  ^  exp[iK.d]  j  «  exp[-iQt]exp[K.Z' (d,t) .K] 

d 

=  exp[-2W(K) ]  'Y  [K..  .K]..  .$n(K,fl)..  .[K..  . K]/n !  (16c) 

n=0 
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Equations  (14a)  and  (16c)  define  ($n(K,fl)},  which  are  designated  as  the  n-phonon 
structure  tensors  and  exp[-2W(K) ] (=exp[-K.W.K] )  is  the  Debye-Waller  factor.  (We 
can  also  relate  the  (In(K,fi)>  to  the  ($n(K,Q)}  by  expanding  exp[-2W(K)],  forming 
the  product,  and  identifying  the  tensor  coefficients.) 

THE  LOW  ORDER  STRUCTURE  TENSORS 

Since  the  structure  tensors  play  a  central  role  in  diffraction  models,  it 
is  worthwhile  to  display  the  expressions  for  the  first  few  in  some  detail.  The 
expressions  are  deduced  by  examination  of  expansions  in  Eq.  (16). 

The  Zero  Rank  Structure  Tensor  In(K>fl)  and  0-Phonon  Structure  Tensor  $n(K,fl) 

*o(K,fl)  =  I0(K,fl)  *  N  J  exp[iK.d]  j  ~  exp[-iflt]  =  a(K)6(fl)  (17) 

d 

defines  the  zero  rank  structure  tensor  I0(K,n)  and  zero  rank  0-phonon  structure 
tensor  $0(K,O) . 

The  Second  Rank  Structure  Tensor  I^(K,0)  and  1-Phonon  Structure  Tensor  $i(K,fl) 

IjfK.fl)  =  N  ^  exp[iK.d]  j  ~  exp[-ifit]  Z(d,t) 
d 

=  N  2  exp [ i K . d ]  j  exp(-iflt] 
d 

2  wq  t-<2nq+l>  +  <nq+l>exp(i0d;t)  +  <nq>exp(-i0c|/t)] 

q 

=  -a(K)<5(fl)W  +  (K,fl)  =  I0(K,ft)(-W)  +  4>1(K,D)  (18a) 

where 

$l(K,D)  =  2  «q  [<nq+l>a(K+q)6(n-Jlq)  +  <nq>a(K-q)5(D+nq) ]  (18b) 

q 
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is  the  1-phonon  scattering  structure  tensor,  Eq.  (14b)  defines  the  (second  rank) 
Oebye-Waller  tensor  W,  and  we  have  used 

a(K±q)<5(-ft+flq)  =  ^  exp[iK.d]  j  ~  exp[-iS2t]exp(±i0d(  t)  (18c) 

d 

The  Fourth  Rank  Structure  Tensor  I<?(K,Q)  and  2-Phonon  Structure  Tensor  $?(K,fl) 

1 2 ( K > ) i , i 1 , h 1 , h  =  N  ^  exp[iK.d]  |  ~  exp[-iflt]Z(d, t) >h  Z(d,t).j>#h> 

d 

=  N  ^  exp[iK.d]  |  expt-iQt] 
d 

/  fwq' 3 i ' , h ' (<2nq'+l>  -  <nq<+l>exp( ifl 'd( t)  -  <nq« >exp(-i0’di t) ] 

q' 

^  twq] i , h  [<2nq+l>  -  <nq+l>exp( i0dt t)  -  <nq>exp(-i0dft)j 

q 

=  a(K)6(0)  [W2 ] i ,  i ' , h 1  , h  "  [W]i,h  [*l(K,fl)  U  •  ,h' 

-  [<&i(K,n)]i(h  CW]  i  * ,  h '  +  ]  i  #  i  * ;  h  ’  ,h  (19a) 

where 

<J>2(K,fl)  =  ^  rq,q'  <<nq+l><nq<+l>a(K+q+q'  )5(0-nq-fiq> ) 

q,q' 

+<nq+lxnq>>a(K+q-q'  )6(ft-Qq+Qq< ) 

+<nqxnq>  +  l>a(K-q+q'  )5(n+gq-f2qi ) 

+<nqxnq'>a(K-q-q'  )5( fl+flq+nq> ))  (19b) 

is  designated  as  the  2-phonon  scattering  structure  tensor. 
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(19c) 


Tq^q*  —  Sq 1  Gq6qGq  ■  / [  ( 2NM )  ^  Qql2q*J 
and 

[w2]i,i\h\h  =  wi,hwi’h’  (19d) 

where  W-j^  is  the  i,h  component  of  the  Debye-Waller  tensor.  (One  might  designate 
(W2 ] i , i ‘ , h ’ , h  as  ^e  fourth  rank  Debye-Waller  tensor.) 

EXPERIMENTAL  EVALUATION  OF  S(K,Q),  In(K,fl),  W,  AND  #n(K,a) 

As  is  well  known,  S(K,fl)  may  be  directly  determined  by  neutron  scattering 
experiments.  It  is  also  clear  that,  in  principle,  one  may  deduce  the  tensor 
In(K,Sl)  from  analysis  of  the  expression  for  S(K,fl)  contained  in  Eq.  (16b). 
Similarly  one  could  deduce  W  and  ($n(K,n))  from  analysis  of  the  expression  for 
S(K ,fl)  contained  in  Eq.  (16c)  or  from  (In(K,fl)>.  Of  course,  in  practice  the 
experimental  determination  of  S(K,n)  is  exacting  and  the  experimental  evaluation 
of  In(K,$2)  or  <5n(K,fl)  for  n>l  may  be  impractical.  Nevertheless,  it  is  important 
to  note  that  the  same  tensors  In(K,Sl)  or  W  and  $n(K ,fl)  appear  in  the  BW 
expression  and  in  the  PW  expression  and  that,  in  principle,  one  can  deduce  the 
tensors  appropriate  for  treating  the  BW  case  from  the  analysis  of  PW  scattering 
experiments. 

THE  STANDARD  EXPRESSIONS 
PW  Cases 

The  standard  form:  electron  scattering  in  liquid  and  amorphous  metals  and 
neutron  scattering 

When  PWs  are  appropriate  basis  functions,  Eqs.  (15a)  and  (15b)  are 
appropriate  and  the  two  forms,  Eqs.  (16b)  and  (16c),  define  the  appropriate 
structure  factor,  i.e.,  the  Van  Hove  dynamical  structure  factor  (ref  1). 
Furthermore,  when  a  PW  basis  is  employed  vanishing  unperturbed  potential, 

(U(x)  =  0)  is  appropriate,  and  thus, 
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W0{k,k')  =  V0 ( k , k ' )  ->  V(K) 


(20a) 


Also,  recalling  Eq.  (13) 

Vn(k,k' )  =  [AA. . .A]V(k,k' )  — >  [KK...K]V(K)  =  Vn(K)  (20b) 

so  that,  employing  Eqs.  (20)  and  (16b),  Eq.  (12)  becomes 

X(K,fl)  =  X(k,k\fl) 

=  ^  Vn(k,k* )*. .  .In(K,fl)..  Vn(k,k')/n!  +  |  W0(k,k’ )|  2  Io(K,fl) 

n=l 


=  |  V ( K ) j  2  {  ^  [KK...K]..  , In(K,fl) . .  .[KK...K]/n!  +  I0(K,fl)> 

n=l 

=  !  V(K)|  *  S(K,fl)  (21) 

Equation  (21)  is  immediately  recognized  as  the  standard  form  for  electron  scat¬ 
tering  in  liquid  and  amorphous  metals  and  neutron  scattering.  That  is,  as  they 
must,  the  BW  expressions  reduce  to  the  standard  forms  in  the  special  case  of 
PWs. 

Approximate  forms 

The  approximate  forms  are  based  on  Eq.  (16c)  for  PWs. 

Markowitz's  approximation  -  Markowitz  (ref  5)  observed  that  the  temperature 
variation  of  the  electrical  resistivity  of  relatively  high  resistivity  amorphous 
metals  could  be  reasonably  represented  by  neglecting  all  but  the  elastic  scat¬ 
tering  term  in  the  expression  for  S(K,fl),  i.e., 

S(K,fl)  a  exp[-2W(K) ]  4>0(K,fl)  (22) 

where  <{>o(K,fi)  is  defined  in  Eq.  (17).  (Generalization  for  anisotropic  materials 
is  obvious.)  The  fact  that  such  a  drastic  approximation  yields  reasonable 
agreement  with  experiment  suggests  a  breakdown  of  conventional  theory  of 
electrical  transport  in  the  relatively  high  resistivity  amorphous  metals 
studied. 
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Sham-Ziman  approximation  -  The  Sham-Ziman  (ref  6)  approach  to  approximating 


the  multiphonon  contributions  to  the  scattering  rate  is  the  most  popular. 
Reference  6  presents  arguments  to  support  this  procedure,  but  the  simolicity  of 
the  resulting  forms  and  their  relatively  good  performance  is  at  the  core  of 
their  wide  acceptance.  One  assumes  that  the  effect  of  the  multiphonon  series  is 
to  cancel  the  effect  of  the  Debye-Waller  factor  in  the  1-phonon  term  in  Eq. 

(16c),  i.e., 

S(K,fl)  ~  exp[-2W(K)  ]  ‘i>0(K,n)  +  K.*-|(K,fl).K  (23a) 

where  <{>g{K,ft)  and  <f>i(K,tt)  are  defined  in  Eqs.  (17)  and  (13b).  For  isotropic 
materials  (e.g.,  amorphous  metals),  Eq.  (23a)  becomes 

S(K,0)  *  exp[-2W(K)  ]  4>0(K,Q)  +  <5>i(K,n)K2  (23b) 

Other  approximations  -  Other  approximations  are  employed.  Sometimes  the 
multiphonon  contributions  are  simply  ignored.  Meisel  and  Cote  (ref  7)  have  com¬ 
pared  an  approximation  proposed  by  Hernandez-Calderon  et  al.  (ref  8)  with  the 
Sham-Ziman  approximation  (ref  6).  Little  basis  for  selection  of  one  over  the 
other  could  be  found  in  fitting  experimental  data. 

Baym's  Form:  The  Nearly  Free  Electron  (NFE)  Case 

Here  we  consider  an  important  special  case  of  PW  scattering.  Strictly 
speaking  PWs  are  not  appropriate  electron  basis  functions  for  the  treatment  of 
crystalline  materials.  Nevertheless,  PWs  are  frequently  a  reasonable  approxima¬ 
tion.  We  shall  refer  to  such  instances  as  NFE  cases. 

The  NFE  procedure  for  treating  scattering  of  electrons  due  to  lattice 
vibrations  in  perfectly  crystalline  solids  is  to  use  PWs  as  approximate  basis 
functions  of  the  unperturbed  Hamiltonian  Hg,  which  has  a  potential  of  the  form, 

^U(x-l)  (24) 

1 
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where  V(x)  =  U(x).  Thus,  W(K)  =  0,  and  Eq.  (13)  becomes 


X(k,k',fl)  =  2  Vn(k,k')*., 

.  .In(K,fi) . .  .Vn(k,k')/n! 

(25a) 

n=l 

=  1  V(K)j  2  {  ^  [KK...K].. 

.In(K,fi)..  . [KK. . .K]/n! 

(25b) 

n=l 

=  i  V(K)|  2 

s’ (K,n) 

(25c) 

where  one  uses 

Vn(k,k')  ~  Vn(K) 

=  [K..  . K] V ( K) 

(25d) 

which  holds  for  PWs.  Note  that  the  elastic  term  (n  =  0)  is  absent  from  Eq. 

(25).  Equation  (25)  is  frequently  applied  in  crystalline  metals  (e.g., 

Reference  2);  we  refer  to  this  form  as  Baym's  approximation  or  the  NFE  approxi¬ 
mation. 

Bloch  Wave  Scattering  in  Perfect  Crystals 
The  standard  form 

The  procedure  for  treating  scattering  of  electrons  due  to  lattice  vibra¬ 
tions  in  perfectly  crystalline  solids  (for  which  an  NFE  model  is  inadequate)  is 
to  use  (approximate)  BW  solutions  of  the  unperturbed  Hamiltonian  H0  as  basis 
functions,  where  H0  has  a  potential  of  the  form, 

^U(x-l)  (26) 

1 

where  V(x)  =  U(x).  Thus,  W(K)  =  0,  and  Eq.  (13)  becomes 

X(k,k’,n)  =  ^  Vn(k,k' )*. .  .In(K,n)..  .Vn(k,k')/n!  (27) 

n=l 

Note  that  Eqs.  (26)  and  (27)  are  identical  to  Eqs.  (24)  and  (25a);  however, 
since  Eq.  (25d)  does  not  apply  to  BWs,  the  forms  in  Eqs.  (25b)  and  (25c)  do  not 
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apply  to  BWs.  As  in  Baym's  approximation,  the  elastic  term  (n  =  0)  is  absent 
from  Eq.  (27);  i.e.,  the  "natural"  Bloch  electrons  are  not  elastically  scattered 
in  perfect  crystals.  Equation  (23)  is  the  appropriate  form  to  substitute  for 
a(K)  in  the  expressions  for  In(K,Q)  in  the  perfect  crystal  case. 

Approximate  forms:  Sham-Ziman  approximation 

The  Sham-Ziman  approximation  (ref  6)  was  originally  introduced  in  the  con¬ 
text  of  perfect  crystals.  *  The  Sham-Ziman  approximation  for  BW  scattering  yields 

X(k,k*,fl)  *  v1(k,k,)*.$1(K,n).v1(k,k’)  (28) 

Thus,  substituting  Eq.  (23)  for  a(K),  one  finds 

( 2rr ) 3  \  I  ea.Vi  (k,k'  )|  2 

X(k,k',fl)  =  2  {<nq+1>6(K+G+q)5(n"V 

c  q 

G,q 

-<nq><5(K+G-q)<5(R+flq)}  (29a) 

=  ((27I)VVC]  |  g(k,k')|  2  {<nq+l>S(K+G+q)6(ft-Qq) 

G  c, 

-<nq>5(K+G-q)5(Sl+r.q) )  (29b) 

Equation  (29)  is  in  the  form  usually  employed  to  treat  electron-phonon  interac¬ 
tion  controlled  scattering  of  Bloch  electrons  in  crystalline  metals.  That  is, 
the  typical  study  of  electron-phonon  interaction  centre1 led  processes  employing 
BW  electrons  implicitly  assumes  perfect  crystallinity  and  invokes  the  Sha~  Ziman 
approximation.  For  example,  Eq.  (29)  gives  the  "first  principles"  form  used  in 
the  work  of  Butler,  Pinski,  and  Allen  (ref  9).  (One  assumes  that  in  Reference 
9,  6{K~q)  stands  for  a  summation  over  RLV  G  on  5(K+G-q).)  Thus,  the  diffraction 
model  origin  and  the  nature  of  the  multiphonon  scattering  approximation  implic¬ 
itly  adopted  in  standard  treatments  of  electron-phonon  interaction  controlled 
processes  is  revealed. 
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ELECTRON- PHONON  SPECTRAL  FUNCTIONS  AND  STATIC  STRUCTURE  FACTORS 


When  dealing  with  electron-phonon  interaction  controlled  phenomena,  it  is 
customary  to  define  various  electron-phonon  spectral  functions,  denoted 
a2(Sl)F(R).  For  example,  the  Eliashberg  function  (ref  10)  is  an  electron-phonon 
spectral  function.  The  a2(0)F (0)  may  be  defined  by  integrating  appropriately 
weighted  diffraction  model  scattering  rates  over  K.  These  electron-phonon 
spectral  functions  are  of  theoretical  importance;  however,  the  procedures 
followed  in  applying  them  often  obscure  the  diffraction  model  origins  of  the 
analyses. 

On  the  other  hand,  when  dealing  with  electron-phonon  interaction  controlled 
phenomena,  it  is  also  possible  to  define  static  structure  factors,  which  would 
be  denoted  S(K),  Im(K),  or  $m(K).  The  various  static  structure  factors  (or  ten¬ 
sors)  may  be  defined  by  integrating  appropriately  weighted  diffraction  model 
scattering  rates  over  n.  The  x-ray  and  resistivity  static  structure  factors  for 
amorphous  metals,  as  discussed  by  Meisel  and  Cote  (ref  11),  are  examples.  The 
static  structure  factors,  which  are  wholly  determined  by  the  lattice  dynamics, 
are  also  of  theoretical  interest  and  emphasize  the  diffraction  model  origins  of 
the  analysis  being  performed. 

CONCLUDING  REMARKS 

A  generalized  diffraction  model  appropriate  for  EW  electrons  as  well  as  PW 
scattering  from  a  dynamic  lattice  has  been  presented.  It  is  shown  that  the 
resulting  model  reduces  to  the  standard  diffraction  model  in  the  case  of  PW 
scattering.  These  findings  unify  the  treatment  of  diverse  phenomena.  For 
example,  neutron  scattering  in  amorphous  solids  and  electron-phonon  interaction 
contributions  to  electrical  transport  in  perfect  crystals  are  conveniently 
treated  in  the  context  of  the  generalized  diffraction  model. 
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Furthermore,  it  has  been  explicitly  demonstrated  that  many  important  stand¬ 
ard  expressions  appearing  in  the  literature  are  obtained  as  approximations  to 
the  diffraction  model.  Thus,  the  diffraction  model  origins  of  the  standard 
expressions  are  revealed  and  the  nature  of  the  implicit  approximations  in  the 
standard  forms  is  clarified. 
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WASHINGTON,  D.C.  20375 


NO.  OF 

NO.  OF 

COPIES 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 

COPIES 

1 

ATTN:  AFATL/MN 

EGLIN  AFB,  FL  32542-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNF 

1 

1 

EGLIN  AFB,  FL  32542-5434 

1 

27709-2211 

METALS  AND  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  DIVISION 

505  KING  AVENUE 

DIVISION  1 

1 

COLUMBUS,  OH  43201-2633 

1 

NOTE: 


PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL, 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES, 


